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We calculate the transverse momentum dependence of the medium-induced gluon 
energy distribution radiated off massive quarks in spatially extended QCD matter. 
In the absence of a medium, the distribution shows a characteristic mass-dependent 
depletion of the gluon radiation for angles 9 < m/E, the so-called dead cone effect. 
' Medium-modifications of this spectrum are calculated as a function of quark mass 

m, initial quark energy E, in-medium pathlength and density. Generically, medium- 
induced gluon radiation is found to fill the dead cone, but it is reduced at large 
£NJ | gluon energies compared to the radiation off light quarks. We quantify the result- 

ing mass-dependence for momentum-averaged quantities (gluon energy distribution 
and average parton energy loss), compare it to simple approximation schemes and 
discuss its observable consequences for nucleus- nucleus collisions at RHIC and LHC. 
In particular, our analysis does not favor the complete disappearance of energy loss 
ji^. effects from leading open charm spectra at RHIC. 
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1. Introduction 



Hadronic jets accompanying heavy quarks c, b differ from light quark and gluon 
initiated jets. These differences can be attributed to the suppression of gluon 
bremsstrahlung from massive charges at small angles 9 < m/E, the dead 
cone effect [1]. Observable consequences of this mass- dependence of the par- 
tonic fragmentation pattern include the softening of the light hadron spectra 
accompanying heavy quark jets [2], and the significant hardening of the lead- 
ing charmed [3] or beauty [4] hadron. Mass-dependent dead cone conditions 
are implemented in the modified leading logarithmic approximation which 
accounts for jet multiplicity distributions [5]. They are also implemented in 
modern Monte Carlo simulations [6] which provide a probabilistic implemen- 
tation of the perturbative part of the parton fragmentation process in the 
vacuum. 



How is this parton fragmentation modified if the produced high-energy parent 
quark propagates through dense QCD matter, as is the case in ultrarelativistic 
nucleus-nucleus collisions at RHIC and at the LHC? As a first step towards 
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addressing this question [7] , several groups [8-12] calculated to leading order in 
energy and for an arbitrary number of medium-induced momentum exchanges 
the modifications to the q — > qg splitting process for energetic light quarks. 
These calculations indicate that medium effects can result in a significant ad- 
ditional energy degradation of the leading hadron which grows approximately 
linear with the density of the medium and approximately quadratic with the 
in-medium pathlength. Recent measurements [13-18] of high-p^ hadroproduc- 
tion and its centrality dependence in Au+Au collisions at ^/snn = 200 GeV 
provide the first evidence [19] for this medium-induced parton energy loss. 
For light quarks and gluons, the formalism was also extended to the angular 
dependence of the medium- modified gluon radiation [20-23]. This allows to 
discuss medium modifications of jet shape and jet multiplicity observables [24]. 

For massive quarks, the formalism of medium-induced gluon radiation is not 
developed to the same extent. Dokshitzer and Kharzeev [25] suggested that 
the dead cone effect also reduces the medium-induced gluon radiation, thus 
resulting in a smaller suppression of leading charmed and beauty hadrons. 
They estimated this effect by multiplying the medium-induced gluon spec- 
trum for massless quarks with a transverse momentum averaged suppression 
factor given in Eq. (4.16) below. However, the combination of vacuum-induced 
and medium-induced gluon energy distributions is known to differ significantly 
from a simple superposition [20] due to the non-trivial interplay of interference 
effects and elastic scattering. Hence, it is conceivable that the medium-induced 
part of the gluon radiation differs significantly from this averaged dead cone 
approximation. Going beyond this approximation may also be needed to dis- 
tinguish mass- dependent final state effects from non-linear modifications of the 
initial gluon distribution for which open charm production may be a sensitive 
probe [26-29]. 

This motivates to parallel for the massive case the calculations of medium- 
induced gluon radiation which exist for the massless case. To this end, two 
groups [30-33] presented detailed calculations of parton energy loss for mas- 
sive quarks. However, these calculations are limited to the average energy loss 
of massive quarks only. Here, we go beyond these studies i) by providing the 
first analysis of the double differential medium-induced gluon distribution as 
a function of transverse momentum and gluon energy and ii) by comparing 
this gluon distribution to the massless limit for a wide parameter range in 
quark mass, in-medium pathlength and medium density. In Section 2 we set 
up the path-integral formalism for parton energy loss. In Sections 3 and 4, 
we summarize the result of our numerical calculations and we provide simple 
analytical arguments for how the mass-dependence shows up in the medium- 
induced gluon radiation. The main results and their relevance for heavy quark 
production in nucleus-nucleus collisions at RHIC and LHC are discussed in 
the Conclusions. 
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2. Medium-induced Gluon Radiation off Massive Quarks: Formalism 



The gluon distribution radiated from a massive quark traversing spatially 
extended QCD matter can be written as [10,20] 
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(2.1) 



Here, u and kj_ denote the energy and transverse momentum of the emit- 
ted gluon, respectively. The Casimir factor Cf — 3 determines the coupling 
strength of this gluon to the massive quark. For numerical results, we fix the 
coupling constant to a s — 1/3 unless stated otherwise. Eq. (2.1) resums the 
effects of arbitrary many medium-induced scatterings to leading order in 1/E. 
Properties of the medium enter (2.1) via the product of the time-dependent 
density n(£) of scattering centers times the strength of a single elastic scat- 
tering er(r). This dipole cross section <r(r) is given in terms of the elastic 
high-energy cross section |a(q)| 2 of a single scatterer in the colour octet rep- 
resentation, 

ff(r) = 2 /(W |a(q)|2 ( 1 - e ^' r )- (2 ' 2) 



A detailed discussion of (2.1) including the physical interpretation of the in- 
ternal integration variables (yi, yi, y, u, £) can be found in Ref. [10]. 

The only mass-dependence of the gluon distribution (2.1) comes from the 
phase factor exp [iq(yi — yi)}, where q is defined as the difference between the 
total three momentum of the initial quark (pi), and the final quark (p 2 ) and 
gluon (k), 

1 1 

x m uj 

q = Pi-p2~k~— — , x = — . (2.3) 



For the abelian case, the same phase factor is known to give the mass- 
dependence of medium- induced photon radiation to leading order in x 1, 
see Ref. [34,35]. Paralleling the derivation of Ref. [10] for massive quarks, we 
have checked explicitly that this phase is the only mass dependence of the 
gluon distribution (2.1). 

In the absence of medium effects, the gluon energy distribution (2.1) reduces 
to the unperturbed radiation I vac associated to the production of a massive 
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quark in the vacuum. We write the full radiation spectrum as the sum of this 
vacuum component and its medium-modification J med ; 

dl dl vac dl mcd , . 

oj = 00 hw . (2.4) 

duj dk duj dk dw dk 



By construction, both the full gluon distribution oj du , as well as the vacuum 

component oj have to be positive for all values of k± and oj. In contrast, 

the medium-induced modification oj can be negative in some part of 

phase space: negative values correspond to a medium-induced depletion of the 
vacuum component. 

In what follows, we calculate the medium-induced gluon distribution (2.1) in 
two different approximations. We limit the discussion to the case of a static 
medium with in-medium pathlength L for which 

n(0 =no6(L-0- (2.5) 



From the analysis of Ref . [36] , we expect that the case of an expanding medium 
can be reformulated in terms of a static medium (2.5) with suitably adjusted 
density n . 



3. Opacity expansion 



The opacity expansion of the gluon distribution (2.1) amounts to an expansion 
of the integrand of (2.1) in powers of [n(£) a(r)] N . Technical details of this 
expansion can be found in Appendix A and in Refs. [10, 11]. 



A. N = vacuum term: the dead cone effect 



In the absence of medium effects, n(£) = 0, the gluon distribution (2.1) reduces 
to the zeroth order in opacity, 

u dl™ _ JI(N = 0) = OsCf h 
duj dk 1 doj dk 1 n 2 



where H(k±) denotes the radiation term associated to the hard parton pro- 
duction, 
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By construction, this is the vacuum term in (2.4). It shows the dead cone 
effect: gluon radiation is suppressed for gluons which are emitted under small 
angles 



m 2 



(3.3) 



The vacuum term (3.2) for the massive case vanishes for k± — > 0, while the 
corresponding massless limit formally diverges like -rr. 

B. N = 1 Opacity correction to the dead cone effect 

We consider a medium of spatially extended QCD matter, modeled as a col- 
lection of colored Yukawa-type scattering centers (A. 6) with Debye screening 
mass \i. To first order N — 1 in opacity, the average momentum transfer from 
the medium to the hard quark is fi. 

Qualitative arguments: Consider a massless quark first. A gluon of energy uj 
decoheres from the wave function of this quark if its typical formation time 
^coh = ~j is smaller than the typical distance L between the production point 
of the parton and the position of the scatterer. Hence, gluon radiation occurs 
if the phase 



To discuss the transverse momentum dependence of gluon emission, we con- 
sider the corresponding kj_-dependent phase accumulated due to scattering of 
the gluon, 



t co h U 



(3.4) 



exceeds unity. 



Here, the characteristic gluon energy is 




(3.5) 
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(3.6) 



where we use the rescaled transverse momentum 



R 2 = 



/i 



2 ' 



(3.7) 



For sufficiently large transverse momentum, the medium-induced gluon distri- 
bution will show the characteristic perturbative powerlaw for gluon production 
in a hard process, 

J rmcd 1 

u;^^oc- for R>1. (3.8) 
dui dn z k 



This behavior can be checked in the N — 1 opacity expansion by taking the 
large- K limit of Eq. (3.14) below. For transverse momentum R < 1, the gluon 
distribution levels off to a constant [24] which depends on uj/uj c . According to 
the condition (3.6), gluons can be emitted only for R 2 > This allows to 
estimate from (3.8) the transverse momentum integrated gluon distribution 
in the region uj > uj c where only R > 1 is relevant, 

dIZ-n f 00 dR 2 uj c . . 

tc ,_^ z o ^ — oc — for uj > uj c . 3.9 

dui Jui/uic K UJ 



This large-a; behavior has been established for the N — 1 opacity approxima- 
tion [23,37]. 

To estimate the mass- dependence of the gluon distribution, we first introduce 
the dead cone factor 



where 

- 9 x 2 m 2 1 (m 2 \ R - , 

m =— R=QcL - (3 - u) 



For the vacuum term (3.1), this factor has the property that 

J rvac J rvac 

= (3 - 12) 



If the same estimate holds for the medium-induced part of the gluon radiation, 
then the distribution (3.8) is depleted for R 2 < max 7, M 2 . For the large-cu 
region, where M 2 > 7, the transverse momentum integrated gluon distribution 
can thus be estimated from (3.8), 
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dlZ" rdR 2 (m 2 \-u,l , _ ,/m 2 -V /3 



This distribution drops off faster (oc ^) than the corresponding massless term 
(3.9). We thus expect a mass- dependent depletion of the medium-induced 
gluon radiation at large gluon energy. We emphasize that the dead cone sup- 
pression is a sufficient but not a necessary condition for this large-u; behavior. 
What is needed is only that mass effects regulate the singularity qona scale 
M. Since the phase space for R 2 < M 2 is small, it is not essential whether this 
regulation occurs by complete extinction of the radiation e.g. via F(R, M), or 
whether the dead cone is filled by a finite non-singular spectrum for R 2 < M 2 . 
We shall find that the latter case is realized. 

Quantitative analysis: We have calculated the medium modification of the 
gluon energy distribution (2.4) to first order in opacity, see Appendix A. In 
the massless limit, our result reduces to the N — 1 opacity result given in 
Ref. [10]. For finite in-medium path-length L, it interpolates between the to- 
tally coherent and totally incoherent limiting cases of the gluon radiation spec- 
trum. In particular, the incoherent limit is an independent superposition of 
three distinct radiation terms: the hard radiation term (3.2) shifted to H(k+q) 
due to elastic scattering, a Gunion-Bertsch radiation term associated to gluon 
production due to a single scattering center well separated from the quark pro- 
duction, and a probability conserving third term (see Appendix A). In terms 
of the dimensionless variables introduced above, the medium-modified gluon 
distribution to first order in opacity is 

dI(N=l) a s C Ff „ T , r i Aq 2 + M 2 )-Uirih(q 2 + M 2 ) 
" =^(2^) / d?± 



duo dn 2 7r 



(q 2 + M 2 ) 

f (R 2 + M 2 ) + (R 2 - M 2 ) (R 2 - q 2 ) 

g2 + + ^ + - 2 + - 2 ^2 _ 4 - 2 -2 



The medium-induced gluon energy distribution is obtained by integrating 
(3.14) over transverse momentum up to the kinematic boundary k — to. In 
terms of the integration variable R 2 , this corresponds to 

dI(N = 1) fRM 2 , 2 dUN = 1) 



du 



[ H W , 9 dlUV = 1 . 



where the integration boundary is written in terms of R = uo c L. The integral 
of (3.15) over gluon energy defines the average medium-induced energy loss 
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(A£ ind > = f E du J I(N , 1} . (3.16) 
Jo auj 



In the L — > oo limit, equation (3.16) coincides with the expression given by 
Djordjevic and Gyulassy [32]. At finite in-medium pathlength L, the differ- 
ences between (3.16) and Ref. [32] may be due to the use of a different density 
distribution n(£) of scattering centers. 



C. Numerical results 



Explored parameter space: The double differential gluon distribution (3.14) 
for massless quarks depends on two parameters only, the characteristic gluon 
energy uo c and the typical transverse momentum (Debye screening mass) fi. 
Presenting our results in rescaled variables R 2 = k^//i 2 and 7 = u c /uj, we 
explore the unrestricted parameter range for lj c and /i 2 . For the transverse 
momentum integrated gluon distribution (3.15), the parameter R = u c L en- 
ters via the integration boundary. We motivate our choice of the value of 
R from a model analysis of high-p^ suppressed hadroproduction in central 
Au-Au collisions at RHIC. The order of magnitude of the suppression is in 
qualitative agreement with the parameter choice R = 2000, u c = 67.5 GeV, 
and n L = 1, see Ref. [23]. We thus choose R = 1000 and a significantly 
larger value R = 40000 for numerical calculations. The latter can be viewed 
as an upper estimate for LHC. If the quark mass is finite, the double differ- 
ential gluon distribution also depends on the effective mass (3.11) which is a 
function of the ratio m/E and of R. We prefer to present all results in terms 
of m/E and R although they appear in (3.14) only in one combination. 

Results: We studied the differences between the medium-induced gluon dis- 
tribution (3.14) radiated off massive quarks, and the corresponding double 
differential gluon distribution for massless quarks 

dI(N = 1) a s C F r g 2 -±sin( 7 g 2 ) 



m ^o oL» dn 2 7r Jo ' q A 



k 2 — q 2 ^ 



x 7 1 2 ' l3/2 - ( 3 - 17 ) 

[l + R 2 + q 2 f - 4/t 2 g 2 



Also, we tested the conjecture that the medium-induced distribution can be 
obtained from the massless expression by multiplying with a dead cone factor 
(3.10), 

dI dcad (N = 1) - dI(N = 1) 

u ; — — = F(k, M) hm LO ; — — — . (3.18) 

dudn 2 dudn 2 
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Fig. 1 shows that the transverse momentum dependence of the medium- 
induced gluon distribution deviates qualitatively from the dead cone approxi- 
mation (3.18). At small transverse momentum, the medium- induced radiation 
does not vanish as for the dead cone approximation (3.18). In contrast, it is 
enhanced compared to the massless case. 



(n L)~ 1 wdl(N=1 )/dwd/c 2 (n L)~ 1 ajdl(N= 1 )/dwd/c 2 




10 5 10 4 io 3 io 1 1 10 i_o_ 2 2 10 5 10 4 10 3 10 2 io 1 1 10 1_0 2 

K K 



Fig. 1. The medium- induced gluon energy distribution as a function of the 
transverse momentum R 2 = k^/jU 2 and for different values of 7 = Q c /uj. Different 
curves correspond to the full medium-induced gluon distribution (3.14) for a mass 
to energy ratio 0.03 of the heavy quark (solid line), the massless limit (3.17) of this 
spectrum (dotted line), and its dead cone approximation (3.18) (dashed line). 

For a more detailed discussion of Fig. 1, we first note that the mass to energy 
ratio ^7 and the parameter R enter the double differential gluon distribution 
only via the effective mass parameter M 2 = see (3.11). For the numer- 

ical values in Fig. 1, we find M 2 = 0.45/y 2 . In accordance with the qualitative 
arguments given in Section 3 B, mass-dependent deviations are seen to be- 
come significant for transverse momentum R 2 < M 2 . In particular, for smaller 
gluon energies u>, i.e. for larger values of 7 = ^, the onset of mass- dependent 
deviations is at smaller values of R 2 . On the other hand, for larger gluon en- 
ergy (i.e. for small 7), the transverse momentum distribution of the radiation 
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spectrum can exceed the Debye screening mass significantly. The condition 
k 2 < 1 provides only a rough estimate for the upper limit on medium-induced 
gluon radiation. 



(n L)~ 1 cjdl(N = 1 )/dcjd/c 2 , R=1000 




K K K 




K K K 

Fig. 2. Same as Fig. 1 but for the larger mass to energy ratios ^ = 0.1 and 
§ = 0.3. 

We next comment on the origin of the mass dependence in Fig. 1 which is in 
qualitative disagreement with the dead cone prescription (3.18). For massless 
quarks, medium modifications of the gluon radiation are known to be deter- 
mined by two competing effects [20]: First, additional medium-induced gluon 
radiation increases the gluon distribution. Second, medium-induced elastic 
scattering shifts emitted gluons to larger transverse momentum and thus de- 
pletes the oc p- vacuum distribution at small transverse momentum. [For very 

energetic gluons which are emitted predominantly at small angle, the second 
mechanism dominates and the medium-induced part of the gluon distribution 
is hence negative for small transverse momentum. This is seen in the plot for 
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7 = 1 in Fig. 1.] For massive quarks, the dead cone effect implies that there 
is no vacuum distribution at small angle which can be depleted due to elastic 
scattering. As a consequence, for massive quarks the second mechanism does 
not apply at small R and the gluon radiation is further enhanced. 

(n L)~ 1 ajdl(N = 1 )/dw 
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Fig. 3. The medium-induced gluon energy distribution (3.15) calculated from 
the full expression (3.14) for a massive quark (solid line), from the massless limit 
(3.17) (dotted line), from the dead-cone approximation (3.18) (dashed line) and 
from the average dead cone factor (3.19) multiplied by the massless spectrum 
(dash-dotted line). 



We have varied the mass to energy ratio ^ and the parameter R in our calcu- 
lation over a wide parameter range (0.001 < § < 0.3 and 1000 < R < 40000). 
Quantitatively, the medium- induced radiation varies significantly with these 
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parameters (see the discussion below). Qualitatively, the effect discussed above 
is generic for the entire parameter space: For transverse momentum R < M, 
medium-induced gluon radiation fills the dead cone. For a significant part 
of the parameter space, it is enhanced compared to the massless case, too. 
For very large mass to energy ratios, however, mass effects limit significantly 
the medium-induced gluon radiation and the radiation from massless quarks 
finally dominates for R < M, see Fig. 2. 

The phase space of the dead cone region R 2 < M 2 is small. For massive 
quarks, neither the vacuum contribution nor the medium-induced contribution 
to (2.4) contain collinear singularities which could enhance the importance 
of this phase space region. As a consequence, parton energy loss off hard 
quarks will be dominated mainly by radiation in the region R 2 > M 2 . For 
large R, however, the medium-induced radiation off massive quarks (3.14) is 
suppressed compared to the massless limit if the ^-ratio is sufficiently large, 
see Fig. 2. This feature dominates the transverse momentum integrated gluon 
energy distribution (3.15), see Fig. 3. In particular, the radiation is depleted 
for large gluon energy, since it increases with the effective mass parameter 
M 2 = |^r-R^2- Fig. 3 indicates that the mass-dependent suppression of (3.15) 
is even stronger than predicted by the dead cone approximation (3.18). 

One may ask whether the correct amount of mass-dependent suppression can 
be estimated from a /t-independent suppression factor F av d c (M) multiplying 
the gluon energy distribution for massless quarks. [Here, the subscript "avdc" 
stands for "average dead cone".] Such a factor could be useful for simplified 
calculations in which k^-dependent information is not available. Paralleling an 
estimate given in Ref. [25] for the dipole approximation [see Eq. (4.16) below], 
we have estimated F avdc (M) by evaluating the dead cone factor (3.7) for the 
characteristic angle under which medium-induced gluon radiation occurs, 9 2 = 



Substituting this average for the dead-cone factor in (3.18) and calculating 
the transverse momentum integrated spectrum, one tends to overestimate the 
mass-dependent reduction of the medium-induced radiation, see Fig. 3. 

Fig. 4 shows the ratio of the medium-induced average parton energy loss (3.16) 
for massive quarks, normalized to the same quantity calculated for massless 
quarks. In general, a finite quark mass is found to reduce the parton energy 
loss and this reduction increases strongly with ^. 





(3.19) 
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<AE ind > (N = 1) (m/E)/<AE ind > (N = ;) (m/E = Q) <AE ind > ( ^ 1) (m/E)/<AE ind > (N ^ 1) (m/E = 0) 




1 o s 1 o 2 10 ; 1 o 3 1 o 2 1 o 1 

m/E m/E 



Fig. 4. The average medium-induced parton energy loss (3.16) for massive 
quarks, normalized to the massless limit, for different values of the density parame- 
ter R = uj c L. Curves are calculated for the full medium-induced radiation (3.14) off 
massive quarks (solid lines), the dead cone approximation (3.18) (dashed lines) and 
the corresponding expression with averaged dead cone (3.19) (dash-dotted lines). 

We now discuss in more detail the main theoretical uncertainty entering (3.16) 
and Fig. 4. Since the double differential distribution (3.14) is calculated in the 
eikonal approximation, it can have support for gluon energies u > E. This just 
indicates the limited validity of the eikonal approximation for small energies 
E. This issue is critical for calculations of the average parton energy loss 
(3.16) which can depend strongly on the phase space limit ■=-<■§- imposed 
on the integral over the gluon energy distribution. In Fig. 5, we illustrate 
this point by showing the gluon energy distribution and the average energy 
loss (3.16) for the parameter set no L = 4, /j, = 500 MeV, L = 4 fm and 
a s = 0.3. This parameter set was used to reproduce the suppression of high- 
p± hadropro duct ion observed in Au+Au collisions at RHIC [38], and it is 
used in the numerical calculation of Ref. [32] together with a charm quark 
mass of m — 1.5 GeV. For the standard upper integration bound oo < E, the 
massive case agrees qualitatively in magnitude and energy dependence with 
the calculation shown in Fig. 1 of Ref. [32] . Small quantitative differences are 
due to the different finite L dependence of the density distributions, as argued 
following Eq. (3.16). However, as seen in Fig. 5, the average energy loss (3.16) 
off massless quarks turns out to be smaller since a larger part of the gluon 
energy distribution lies above the kinematic cut. We have made no attempt 
to remedy this possibly unphysical behavior. We simply emphasize that any a 
posteriori modification of the large-u; tail of ojj- entails significant theoretical 
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uncertainties. It is an open problem of obvious importance to include finite 
energy constraints on the level of the radiation spectrum (i.e. on the level 
of multiple scattering Feynman diagrams) rather than on the level of the 
integrated parton energy loss. 



wdl(N= 1 )/dw 
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Fig. 5. For finite quark energy E, the normalized average energy loss (rhs) de- 
pends significantly on the kinematic boundary up to which the gluon energy distri- 
bution (lhs) is integrated. This entails significant uncertainties, which are discussed 
in the text. Parameter values are taken from Ref. [32]. 



4. Dipole approximation (multiple soft scattering) 



If the medium provides a large number of soft momentum transfers, rather 
than a few harder ones, then the projectile performs a Brownian motion in 
transverse momentum. This case can be studied in the saddle point approxi- 
mation of the path-integral (2.1), using [34,39] 

n(0 <r(r) ~ \ q(£) r 2 . (4.1) 



Here, q(£) is the transport coefficient [8] which characterizes the medium- 
induced transverse momentum squared {qj_ ) transferred to the projectile per 
unit path length A. In the approximation (4.1), the path integral in (2.1) 
is equivalent to that of a harmonic oscillator which allows for an explicit 
calculation. 
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A. The Vacuum Term shows the Dead Cone Effect 



For the transport coefficient q(£) = qO(L — £) of a static medium of length L, 
we evaluate the gluon distribution (2.1) by splitting the longitudinal integrals 
into three parts [35], 

dl dh dh die 

du dn 2 du dn 2 duj dn 2 du dn 2 

rh rL rL /*oo /*oo /*oo 

= / dyi dy t ...+ dyi dy t . . . + dyi dyi... (4.2) 

JO Jyi JO JL JL Jyi 



In complete analogy to the calculation for the massless case [20], one can 
shown that the term J 6 does not depend on the medium and takes the form 

dh = as C F kj (4 3) 

^dud^ ir 2 (k 2 ± + x 2 m 2 ) 2 ' 1 ' ' 



This expression coincides with the vacuum term (3.1). In the multiple soft 
scattering separation, there is hence a simple separation of the gluon distribu- 
tion (2.4) into the vacuum term (4.3) and the medium-induced contribution 
I 4 + J 5 . The latter vanish in the absence of a medium, q — 0. 



B. Medium-induced gluon radiation 



We now discuss the medium-induced part of the gluon distribution (2.4) which 
in the multiple soft scattering (mss) approximation takes the form 

_ dl mss _ _ dh | „ dh ^ 
du dn 2 du dn 2 du dn 2 



The terms J 4 and I 5 are defined in (4.2) and are given explicitly in Appendix B. 

Qualitative arguments: Consider first the qualitative behavior of (4.4) for the 
case of a massless quark. The energy and transverse momentum scales which 
determine medium-induced gluon radiation can be estimated in analogy to the 
discussion given for the N — 1 opacity case. The phase accumulated by the 
gluon due to multiple scattering is given in terms of the characteristic gluon 
energy u c [40], 

/ k 2 \ q L u c 1 . T o , , _ N 
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This phase should be larger than 7 ~ 0(1) for the gluon to decohere from the 
massless quark. Here, we have used that the gluon carries a typical transverse 
momentum squared of order (kj_) ~ q L. We express the transverse momentum 
kj_ in units of this characteristic scale, 

„' = jj. (4.6) 



For the case of many soft scattering centers, several centers act coherently if the 
formation time of the gluon exceeds its mean free path. The typical transverse 
momentum accumulated by the gluon is then ~ (q\) ^j™- = qt {orm while its 
formation time is tf orm ~ Medium-induced radiation thus relates transverse 
momentum and gluon energy, 



c 2 ~ 



'qoj or in dimensionless units: k —\- — • (4.7) 




On the other hand, the number of scattering centers which act coherently 
is N coh = ^ nL . For iV coh > 1, the medium-induced gluon radiation in the 
multiple soft scattering limit reads 

J rmss T Jrlscatt 

u a m =° ~ - u m=0 , (4.8) 
duj dn 2 A N coh duj dn 2 ' 



where 1^=0* is the radiation off a single scattering center and the total number 
L/Xoi scattering centers along the path of extension L is reduced by coherence 
effects to the effective number 

W * ~ - - (4.9) 



AiV coh V U) K 2 

If the gluon energy is small, uj < cu c , then it follows from (4.7) that k 2 < 1. 

^ylscatt 

The transverse momentum integral over u d ™^ 2 gives an energy-independent 
function and we find from (4.8) and (4.9) that 

A rmss J r—r- 

W %2oc-^ for u<u c . (4.10) 

duj \N coh V uj 

^ylscatt 

For k > 1, on the other hand, the spectrum characterized again 
by the perturbative tail (3.8). Combining this information with Eqs. (4.8) and 
(4.9), one obtains 

uj , ~ / — r ~ — for uj > uj c . (4.11) 

duj Juj/uc K b \UJ J 
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Both limiting cases, (4.10) and (4.11) agree with the results of the full calcu- 
lation [23] in the multiple soft scattering limit for massless quarks. 



To estimate the mass-dependence of the gluon distribution, we parallel the 
discussion of Section 3 B. We introduce the dead cone factor 



2 \ 2 

K 



* M > = fe = sra? • (412) 



where 



„, 9 x 2 m 2 l(m 2 \ ^u 2 „ r tl 

" = T = iN ! 5' fl = " cL - (413) 



Restricting the ^-integration of the perturbative tail by this mass term, we 
find 

dl mss foo 1 

u—r / — pT ~ T77 ■> for cj>w c . (4.14) 

duo Jm* k 6 M 4 v y 



In close analogy to the estimate for single hard scattering in Section 3, the 
large-a; region drops off significantly faster than for the massless case (4.11). 
We thus expect a mass-dependent depletion for large uo. 



C. Numerical results 



Explored parameter space: In analogy to the study in Section 3 C, we present 
all results in rescaled variables k 2 and 7. The parameter R = oo c L is explored 
for the values 1000 and 40000. See Section 3 C for further details. 



Results: We have calculated the transverse momentum dependence of the 
medium-induced gluon distribution (4.4) in the multiple soft scattering ap- 
proximation. We have compared this distribution to the massless limit and to 
the dead cone approximation 

J rmss J rmss 



As seen in Fig. 6, the results are in qualitative agreement with those obtained 
in the opacity expansion (see Fig. 3 and Section 3 C for discussion). In partic- 
ular, we find that the dead cone is filled by medium-induced gluon radiation. 
This feature persists to larger mass to energy ratios (data not shown) in qual- 
itative agreement with the opacity expansion in Fig. 2. 
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Fig. 6. The medium- induced gluon energy distribution as a function of the 

k 2 

transverse momentum k — and for different values of 7 = u> c /u, calculated in 
the multiple soft scattering approximation. Different curves correspond to the full 
medium-induced gluon distribution for a mass to energy ratio 0.03 of the heavy 
quark (solid line), the massless limit of this spectrum (dotted line), and its dead 
cone approximation (4.15) (dashed line). 



The transverse momentum integrated gluon energy distribution calculated 
from (4.4) is shown in Fig. 7. This integral is dominated by the phase space 
region k 2 > M 2 , and thus does not depend strongly on the finite medium- 
induced radiation inside the dead cone region. Compared to the massless case, 
the large-w tail of the gluon energy distribution is depleted with increasing 
m/E-xaXio. The approximation of this effect by the dead cone approximation 
(4.15) tends to underestimate this depletion. These findings are in qualitative 
agreement with the results reported for the opacity expansion in Fig. 3 and 
with the qualitative expectations based on (4.14). 

We have also tested numerically the approximation of Dokshitzer and 
Kharzeev [25]. These authors replaced the transverse momentum dependent 
dead cone factor in (4.15) by an average expression, evaluated at the char- 
acteristic angle Q\ ~ V^ 1 ^- under which medium-induced gluon radiation 



18 



occurs on average, 



DK 



1 m? R 



1 + V2 £ 2 7 3 / 2 



(4.16) 



This approximation allows to mimic the qualitative trend but is found to 
overestimate the depletion of the large-cu tail significantly, see Fig. 7. 



cjd / d co 



: R = 1 O O O 



R = A O O O O 



CO / CO , 




m > C 



Fig. 7. Same as Fig. 3, here calculated in the multiple soft scattering approxi- 
mation. 
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Fig. 8. Same as Fig. 4, here calculated in the multiple soft scattering approxi- 
mation. 



From the gluon energy distribution, we have calculated the average parton 
energy loss according to (3.16). As seen in Fig. 8, a finite quark mass reduces 
parton energy loss significantly for sufficiently large mass to energy ratios, 
m/E > 0.1 say. For smaller mass to energy ratios, there is some parameter 
range where our numerical results indicate the opposite effect. However, similar 
to the case of the opacity expansion, this may be an artifact of the phase space 
constraint in the definition of the average energy loss (3.16) [see the discussion 
of Figs. 4 and Fig. 5 above] . To explore this theoretical uncertainty, we have 
paralleled the logic of Section 3. We have calculated in Fig. 9 the gluon 
energy distribution (4.4) for the set of parameters which reproduce the nuclear 
modification factor for central Au+Au collisions at RHIC [23]. Integrating this 
gluon energy distribution for uo < E, we find again that the average parton 
energy loss for massless quarks can be larger than that for massive ones, simply 
because the kinematic boundary uj < E cuts off the more pronounced large 
cu-tail of the gluon energy distribution for m — 0. For the reasons given in 
Section 3 C, we conclude that this motivates to go beyond the high-energy 
approximation and to include finite energy constraints in the calculation of 
(jJ^t- rather than to impose them a posteriori in the integral over w^. 
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Fig. 9. Same as Fig. 5, here calculated in the multiple soft scattering approx- 
imation for parameter values R = 2000, uj c = 67.5 GeV and m = 1.5 GeV. This 
figure allows to illustrate the uncertainties in calculations of the average parton 
energy loss for finite energy quarks, see text. 



5. Conclusion 



Early studies demonstrated that the transverse momentum spectra of charmed 
and beauty hadrons in nucleus-nucleus collisions depend strongly on the as- 
sumed parton energy loss [41-43]. More recently, the quark mass dependence 
of this effect was argued to give access to the detailed mechanism of medium- 
modified parton fragmentation [25,30-33]. To better assess this mechanism, 
we have presented here the first study of the transverse momentum and energy 
dependence of the medium-induced gluon radiation off massive quarks. 

Our calculation is based on the path integral formalism (2.1) which resums 
medium-modifications to leading order in 1/E and to all orders in opacity. 
We have employed two approximation schemes which model the medium- 
dependence of the hard parton fragmentation as a series of many soft gluon 
exchanges or as a single semi-hard momentum transfer. Despite these different 
physical pictures, both approaches lead to comparable results, in agreement 
with an earlier study of the massless case [23]. 

We find that medium-induced gluon radiation generically fills the dead cone 
9 < m/E [see Figs. 1, 2 and 6]. However, in comparison to the transverse 
momentum integrated gluon energy distribution radiated off massless quarks, 
the radiation is depleted at large gluon energies [see Figs. 3 and 7]. The 
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average parton energy loss results from an interplay of both effects and tends 
to be smaller for massive quarks than for massless ones [see Figs. 4 and 8]. 
However, for sufficiently small parton energies E, the kinematic boundary 
uj < E can lead to the peculiar case that the average parton energy loss is 
larger for massive quarks than for massless ones [see Figs. 5 and 9] . We have 
argued that this behavior may be unphysical and limits the application of the 
formalism of Refs. [8-12] to sufficiently large quark energies, see discussion at 
the end of Section 3 C. 

As of today, the only experimental information about open charm production 
in Au+Au collisions at RHIC is the prompt electron spectrum measured by 
PHENIX at v /s^ = 130 GeV [44]. It was argued that these data do not 
indicate a significant parton energy loss for charm quarks. But it is equally 
true that they do not constrain parton energy loss significantly: First, the 
experimental errors on the prompt electron spectrum are still large. Second, 
for the measured transverse momentum range p^ cctron < 3 GeV, the correla- 
tion between the transverse momentum of the electron and of the charmed 
hadron is very weak. Both complicate any conclusion about the medium- 
induced degradation of p^ iarm . In addition, the moderate values of p c ^ &rxa ac- 
cessed by PHENIX correspond to charm quarks which move rather slowly and 
thus may turn into hadrons prior to leaving the collision region. In this case, 
the energy degradation of charmed hadrons would get contributions from their 
hadronic cross sections and the formalism employed here has to be revisited. 

The experimental information from RHIC is expected to improve soon. First, 
the higher statistics of future runs will allow to measure the prompt electron 
spectrum in a wider p± range. Second, the topological reconstruction of the 
hadronic decay of charmed hadrons (D° — > K~tt + ) should provide a more di- 
rect measurement out to significant transverse momentum. On a longer time 
scale, open charm measurements at LHC will further extend this p± range 
to p®° ~ 15 GeV [45]. Also, the energy loss of beauty quarks is expected 
to become accessible at LHC via high- mass dimuon and secondary J/\P pro- 
duction [46]. Despite the significant uncertainties of our calculation discussed 
above, there is one conclusion which we can draw: parton energy loss is re- 
duced by mass effects, but for realistic parameter values it remains sizable. 
Thus, our study favors a medium- induced enhancement of the D / n° ratio at 
sufficiently large transverse energy but we still expect the nuclear modification 
factor for D° to lie below unity. 

Acknowledgment: We thank R. Baier, A. Dainese, K. Eskola, H. Honkanen, 
A. Morsch, G. Rodrigo and J. Schukraft for helpful discussions. 
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A. Gluon energy distribution to first order in opacity 



In this appendix, we give details for the calculation of the zeroth and first 
order in opacity of the gluon energy distribution (2.1). We start by expanding 
the path integral in Eq. (2.1) in powers of n(£) <x(r), 



1C(r,W,r,yi) = £o(r,yi;r,yi) 



Z Z 

dpi dp 2 /Co(r,j/,;pi,£i) 

1 1 

x 2 a (pi) mpi> Pa, 6) g °"(P2) Mp2, 6; r > j/j) • 



(A.1) 



Here, the free path integral /Co is of Gaussian form 



2iri(yi - y{_ 



exp 



iw r-r 



2 (j/j - y/) 



(A.2) 



We expand the integrand of (2.1) to first order in cr(r). The longitudinal 
integrals in (2.1) are regularized in intermediate steps of the calculation, as 
explained in Ref. [10,35]. We note that the N-th order of (2.1) involves 2A^ + 1 
terms only. This makes it straightforward to obtain explicit expressions to 
high orders in opacity. 

Inserting (A.l) into (2.1), one finds to zeroth order opacity term (3.1). To first 
order in opacity, one finds 



uj- 



dI(N = 1) a s C F 



dcudk 



7T^ 



J dC l L/„M2 



X 



LQi - sin(LQi) 
(27r) 2|alqjl [(k_L + q) 2 + x 2 m 2 f 
(k ± + q) 2 k ± -(k ± + q)" 



(kj_ + q) 2 + x 2 m 2 k 2 ± + x 2 m 2 



(A.3) 



where the transverse energy of the scattered gluon is Q\ = ^ k±+q ^ +:r m . 
Incoherent limit 

In the totally incoherent limit, the first order opacity term (A.3) has a prob- 
abilistic parton cascade interpretation: 



dI(N = 1) 



lim lo 

l^oo duj akj 



noL=const 



7i A 



(no L) 



dq 



|a(q) 
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x [-H(k ± ) + H(k ± + q) + R(k ± , q)] . (A.4) 



Here, H(k±) is multiplied by the probability that the hard parton interacts 
with the medium; the minus sign ensures probability conservation by reducing 
the corresponding weight of the N = vacuum contribution (3.1). For a 
scattering center far away from the point of quark production, interaction 
with the medium gives rise to two processes: i) The vacuum radiation term 
H(k±) is shifted to H(k± + q) due to medium-induced elastic scattering by 
a transverse momentum q. ii) Gluons are produced due to bremsstrahlung off 
the far distant scattering center. This leads to the Gunion-Bertsch radiation 
term for massive quarks, 

*(k„ q) - ^ + (1 » 2 2k ± .(k ± + q) 



[(k_L + q) 2 + x 2 m 2 } 2 [(k_L + q) 2 + x 2 m 2 } [k 2 + x 2 m 2 } 
+ [k 2 +x 2 m 2 ] 2 ' (A ' 5) 



For realistic, finite in-medium pathlength, both effects combine in the specific 
interference pattern (A.3). 

N = 1 opacity term for Yukawa-type scattering potential: 

We have studied (A.3) for arbitrary in-medium pathlength for a Yukawa-type 

elastic scattering center with Debye screening mass /i: 

M q >l 2 = M 2 ^(?W < A - 6) 



Shifting in (A.3) the relative momentum integration q — > q + k and doing the 
angular integrations, one finds 



UJ 



dI(N = 1) _ a s C F n fi 2 
dui dk 2 



roc 

Jo 



.LQi-sin(LQi) q 2 



n uj Jo (.)' q 2 + x 2 m 2 

/i 2 (k 2 + x 2 m 2 ) + (k 2 - x 2 m 2 ) (k 2 - q 2 ) 



x- 



{k 2 + x 2 m 2 ) {m 2 + k 2 + q 2 ) 2 — Ak 2 q 2 



3/2 



(A.7) 



where Q[ = q + 2 ^ m ■ From this expression, we find Eq. (3.14) by rescaling the 
variables uj, k 2 and m 2 to the dimensionless variables 7, R 2 and M 2 defined in 
eqs. (3.4), (3.7) and (3.11) respectively. 
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B. Gluon energy distribution in the dipole approximation 



Here, we give the full expression for the medium-induced part of the double dif- 
ferential gluon distribution (4.4). The calculation is done in complete analogy 
to the calculation of the massless case [20], but keeping the mass-dependence 
of (2.1). In the rescaled dimensionless variables introduced in Section 4, one 
finds 



dh 



a, C 



did dn 2 



-^ 7 2 2Re I' dt / W^^exp 

71 JO Jt 



K 



4(£> 4 - iA 4 B 4 ) 



x 



%Al Ba k 2 



±A 2 D 4 



(D 4 - iA 4 B 4 ) 3 (D 4 



iA 4 B 4 ) 2 



(B.l) 



where 



« = (B.2) 



and 



Aa 



4 7 sin [Q(t - t)} ' 



Ba = COS [Q(t - t)] , Da 



i) . (B.3) 



The term I 5 in Eq. (4.2) takes the form 



dh 



CO 



a s Cf 



did dn 2 



^ 7 2 Re / dte^ M2 ^ - 

71 JO \ K 



— IK 



x-^exp 



2 + M 2 J 

4A; B 5 



(B-4) 



where 



A 5 = - ^t-t , B 5 = cos [<TLt] . (B.5) 
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